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Abstract. Disformal theories of gravity are scalar-tensor theories where the scalar couples 
derivatively to matter via the Jordan frame metric. These models have recently attracted 
interest in the cosmological context since they admit accelerating solutions. We derive the so¬ 
lution for a static isolated mass in generic disformal gravity theories and transform it into the 
parameterised post-Newtonian form. This allows us to investigate constraints placed on such 
theories by local tests of gravity. The tightest constraints come from preferred-frame effects 
due to the motion of the Solar System with respect to the evolving cosmological background 
field. The constraints we obtain improve upon the previous solar system constraints by two 
orders of magnitude, and constrain the scale of the disformal coupling for generic models to 
J\4 > 100 eV. These constraints render all disformal effects irrelevant for cosmology. 
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1 Introduction 

The acceleration of the cosmic expansion is one of the biggest mysteries in modern cosmology 
and theoretical physics. The search for the underlying driving mechanism, dubbed dark 
energy, has previously prompted the study of light cosmological scalar fields as one potential 
candidate (see [1-3] for reviews). Theories where a scalar couples to matter, scalar-tensor 
theories, can be seen as modifications of general relativity (GR); these have been the subject 
of a considerable amount of recent research. Among the plethora of models in the literature, 
those that contain screening mechanisms [4, 5] are particularly well-studied due to their 
ability to decouple solar system scales from cosmological ones. This circumvents the need 
for fine-tuning the model parameters to unnaturally small numbers, as the modifications can 
still be a dominant influence in the cosmological background while satisfying solar system 
bounds. Scalar-tensor theories are most commonly realised through a conformal coupling of 
the scalar cj) to matter via a metric relation of the form 

£lfj,u = A ((f>)g/j,i/, (1-1) 
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where g^ obeys equations of motion derived from the Einstein-Hilbert action (Einstein-frame 
metric), while g^ is the physical metric whose geodesics free-falling matter follows (Jordan- 
frame metric). However, this is not the most general relation one can write down. Bekenstein 
[6] showed that the most general relation that preserves causality is 

9nv = C((f>, X)g^ v + D((f>,X)dn<j>d v (t> where X = —g^d^d^. (1.2) 

The addition of this coupling to any theory does not introduce a ghost by itself. 1 When 
C and D depend on (f> only the theory is a subset of the Horndeski class 2 [7, 8], and when 
D depends on <f> and X the theory fits into the “beyond Horndeski” class [9]. The latter 
case contains hidden constraints that render the equations of motion second-order [10]. The 
term proportional to D is known as the disformal coupling and, in contrast to the first 
term—the conformal coupling—it has only recently begun to be studied in detail [11-26]. 
Until recently, most studies have focused on the cosmology of these theories and, unlike the 
conformal coupling [27-34], little attention has been given to astrophysical and solar system 
effects. Furthermore, Ref. [35] claimed that disformally coupled theories have their own 
screening mechanism. Ref. [21] has examined the local behaviour of these theories and found 
that no such screening mechanism exists. He further found that local objects only source 
scalar gradients (and hence modifications of gravity) when the background cosmological scalar 
(j)Q(t) is accounted for. In this case, the sourcing is proportional to <f >o and <f>o and, while there 
are no non-linear screening mechanisms, the modifications can be screened if the cosmology 
is such that these time-derivatives are small. In this case, the effects on matter are screened 
everywhere, in contrast to non-linear mechanisms that only hide the scalar field locally. Note 
that the natural scale for the background evolution is the Hubble scale, so that at the present 
time we expect <f>o ~ Hofo, <j >q ~ unless the model is fine-tuned. 

In this paper, we rigorously calculate solar system observables for a disformal theory 
with canonical scalar action. In particular, we calculate the parameters appearing in the 
parametrised post-Newtonian (PPN) metric. The PPN formalism [36] is a general framework 
for testing alternate theories of gravity in the solar system. Deviations from GR are encoded 
in 10 parameters (one Newtonian and nine post-Newtonian) that appear multiplying various 
terms in the metric. Solar system experiments constrain these parameters either directly or 
in combinations (see [37] for a review) and so one can constrain a wide class of theories using 
the same experimental data. 

Calculating up to the first post-Newtonian order (2PN), we confirm that there is no 
screening mechanism, and that at all local effects are proportional to cj> q. Thus, the findings 
of Ref. [ 21 ] hold to higher order as well. As a diffeomorphism-invariant theory derived from 
a Lagrangian, disformally coupled scalar-tensor theories are semi-conservative. The non-zero 
PPN parameters are the Eddington light bending parameter 7 , the amount of non-linearity 
in the held equations f3, and the two preferred -frame parameters aq and a-j,. In scalar-tensor 
theories, the evolution of the cosmological background held provides a preferred 4-vector, 
singling out the cosmological (CMB) frame as preferred frame. When the scalar is coupled 
disformally to matter this leads to signihcant non-zero aq, <37 ■ All of the PPN parameters 
are proportional to the same quantity, which is a combination of the time-derivative of the 
cosmological scalar and the disformal coupling scale A. The tightest constraint comes from 

1 This is assuming that the signs of C and D are chosen appropriately. 

J Unless non-Horndeski terms are present in the action S[g t ,^. <p\. 


- 2 - 



ot 21 which is constrained by local tests to be of order 10 -7 or smaller. For generic (non- 
fine-tuned) models, this constrains the coupling scale to A/Hq > 3 X 10 3 . Equivalently, 
normalizing the scalar held to the Planck mass, the coupling is constrained such that AT = 
(MpiA) 1 / 2 > 10 2 eV. This means that any observable consequences of the disformal coupling 
in the cosmos are negligibly small. Moreover, it rules out the disformal coupling as a source 
of the cosmic acceleration. Apart from extreme fine tuning to make 4> o, 4 >o vanishingly 
small today (which will most likely not allow for an accelerating solution without vacuum 
energy), the only way to evade these constraints is to augment the model with the well-known 
chameleon or Vainshtein screening mechanisms. 

The paper is set out as follows: We begin in section 2 by introducing disformal gravity 
theories and give a brief introduction to the PPN formalism in section 3. The calculation of 
the PPN parameters for the minimal disformal model, which is merely chosen for clarity of 
the presentation, is worked out in detail in section 4. In section 5 this is used to place new 
constraints on Hq/A and the first derivative of the disformal factor for a well-studied expo¬ 
nential model. The constraints on the disformal mass-scale AT are two-orders of magnitude 
stronger than previous solar system constraints [21] but the constraints on the first derivative 
of the disformal factor are relatively weak and values as large as 10 3 are not excluded. We 
conclude in section 6. In Appendix A we show that additional contributions to the results 
in the most general disformal theory are sub-leading and hence our results of section 4 and 
5 hold for a much wider class of conformal/disformal models. 


2 Disformal Theories of Gravity 


Canonical disformal gravity theories are described by the scalar-tensor action: 


5 = 


d 4 *^L 


V^9 


Mp\ 2 R(g) i 


- g - VW 


+ 'S'm [g 


Hv J ? 


( 2 . 1 ) 


which describes a massless spin-2 graviton and an additional scalar degree of freedom in the 
Einstein frame. The modifications of general relativity arise due to a coupling of the scalar 
to matter via the Jordan frame metric 


9in> 




g^u + 


b 2 {4>) 

A 2 


dpcpd^cj) 


( 2 . 2 ) 


The dimensionless function A((j>) describes the conformal coupling via the Einstein frame 
metric g while the dimensionless function B(cf>) controls the strength of the derivative 
interaction, which is known as the disformal coupling. Note that we have normalized <f> to 
be dimensionless. In this paper, we will only study the case where the coupling functions 
depend solely on the scalar field, 4>. Our leading constraints will also apply to models that 
generalize this to B((j),X). The mass scale A controls the amplitude of disformal effects, with 
smaller values of A leading to larger effects. 

In what follows we will write M p 2 = (87rG) _1 but it is important to note that G is not 
necessarily equal to Newton’s constant as measured by local experiments and, indeed, we 
will see later that it is not. The equations of motion in the Einstein frame are [21] 


Gfj, v — 8nG (Tm/ii/ + 

1 - 2XB l □</» = - 8vr aGT n 


(2.3) 


A 2 


A 2 
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(2.4) 


o /'"v d2 

—- P+mTiS'd^t + V(<f>), 


where X = —1/2 g^d^cftd^cj) and we have defined 3 

d In A((p) 


a<h{®) = 


d (/> 


and /^(</>) = dln d ^^ ' 


(2.5) 


Here, T£ u = is the energy-momentum tensor of the matter fields defined with 

respect to the Einstein-frame metric and 


T^nv = 


1 


8 vr G 


- 9liV ( + V{<f>) 


( 2 . 6 ) 


is the energy-momentum tensor of the field. Note that these are not separately covariantly 
conserved quantities due to the coupling of the scalar to matter, only their sum is. Equation 
(2.3) is simply Einstein’s equation, which is a consequence of working in the Einstein frame. 
Equation (2.4) is the equation of motion for the scalar and this is where the modifications 
of gravity become apparent. Note that the physical metric is the Jordan frame metric (2.2); 
it is this metric that governs the motion of test particles. The energy-momentum tensor 
of matter appears in equation (2.4) and so the field is sourced by any non-zero matter 
distribution characterised by T mfiu . In non-relativistic systems, the gravitational field is 
sourced mainly by matter and not the scalar, in which case the solution of (2.3) is identical 
to the GR solution. The physical metric, and correspondingly the motion of test particles, 
then deviates from the GR prediction. 

The goal of this paper is to calculate the physical (Jordan frame) metric to post- 
Newtonian order. The post-Newtonian predictions of purely conformal scalar-tensor theories 
have been well-studied (see [37, 38] and reference therein) but, to date, the post-Newtonian 
behaviour of disformal theories has yet to be derived. We therefore focus on the disformal 
part of (2.2), setting A = 1 and a$ = 0. Furthermore, we set B(cj>) to be constant (fy = 0), 
and absorb its value into A to set B = 1. As we will see below, the presence of a disformal 
coupling significantly complicates the calculation compared with the pure conformal case and 
so we will work with this minimal model first in order to make the calculation as simple as 
possible. In fact, in Appendix A we show that the general model yields no new constraints 
compared with the minimal model. This is because the conformal and disformal couplings 
are independently constrained by different PPN parameters, so that any interaction of the 
two effects is highly suppressed. Ref. [21] has examined the Newtonian behaviour of these 
theories and has shown that the disformal terms are only active when one accounts for the 
fact that the space-time is asymptotically FRW and not Minkowski and that the disformal 
terms are sourced by the time-derivative of the cosmological scalar <f>Q. The constraints for 
the non-minimal model are then essentially the same when setting B —> -B(<^o(^o))- 

The minimal model is described by the Jordan frame metric 




df/pd^cj) 

9 ^ H-^ 2 — 


(2.7) 


The inverse Jordan frame metric is given by 


“/XT' 


2X 


= 9^- T A i-TT) 


A 2 


( 2 . 8 ) 


3 We define the derivative of the conformal and disformal factors as to avoid any confusion with the 

PPN parameters cu, /?. The parameter was called 7 in [21, 22]. 
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3 The PPN Formalism 


With the exception of compact objects such as black holes and neutron stars, astrophysical 
objects move with non-relativistic velocities v and one can solve Einstein’s equations by 
expanding in appropriate powers of v/c. In what follows, we will work in units where c = 1 
and describe the expansion in powers of l/c 2m as being O(rri) for brevity. In this sense, the 
solution of the equations to order (v/c) 2 is 0(1), the solution to order (v/c) 5 is 0(1.5) and the 
solution to order (v/c) 4 is 0(2). We will refer to the solutions at 0(1) as 1PN solutions and 
0(2) as 2PN. Velocities v are 0(0.5), while the Newtonian potential is 0(1). Furthermore, 
time-derivatives add a power of v/c, while spatial derivatives do not. The leading effect in 
the disformal model comes from the time derivative of the background field </>o, which we 
take to be 0(0) (in the CMB rest frame). 

The PPN framework parametrises the solution of Einstein’s equations, or more generally 
the Jordan-frame metric, to 2PN as 


500 = -1 + 2U - 2/3U 2 - 2& w + (2 7 + 2 + a 3 + Ci - 2£)$i + 2(3 7 - 2/3 + 1 + C 2 + C)$ 2 
+2(1 + C3)^3 + 2(3 7 + 3(4 — 2£) < I > 4 — (Ci — 2£)-4 — (a\ — a 2 — a 3 )w 2 U — a 2 w l w j Uij 
+(2a 3 - ai)w l Vi, (3.1) 

goi = -^(4 7 + 3 + «1 - a 2 + Cl - 2C)JC - -(1 + a 2 - Cl + 2C)W* - 7 («i - 2a 2 )w L U 


- a 2 w j Uij, 

Qij = (1 + 2 7 P)<%, 


(3.2) 

(3.3) 


where the 1PN potentials are 


U = Gn / d 3 a? 


p(t, x') 


x - x' 


and Uij = Gn / d 6 x 


3 ^,p(t,x')(x - x')i(x - x')j 


\x — x 


'/|3 


(3.4) 


and the 2PN potentials are 


$ lS G N f d 3 S -Kx-> 2 (x') 


x — X’ 


t, 2 = G N [ 




$4 = G n / d 3 x 


\x — X 


t/\ ’ 


\x — X’ 
p(x') 


X — X‘ 


VI ’ 




\x — X ' 




= Cn 


d 3 x / d 3 x „ p(x')p(x")(x~x') ( x' - X” 


\x — X' 


\x — X 
x — x" 


V |3 


\x — x" 


and 




x — x' 


(3.5) 


Note that p(t, x) is the rest-frame mass density of the fluid, which is important when trans¬ 
forming from a given set of coordinates to PPN. The PPN potentials are defined with respect 
to the locally measured Newton’s constant Gn, which is not necessarily equal to the constant 
G appearing in the action. Instead, Gn is defined such that it brings goo at 1PN order into 
the simple form —1 + 2 U. The vector w r , taken to be 0(0.5), is the speed of the PPN coor¬ 
dinate system relative to the mean rest frame of the universe and is taken to be constant in 
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space and time, since it varies over time scales much longer than Solar System time scales. 
The PPN metric is not just a parametrisation of the solution but also a gauge choice: the 
coordinate time is equal to the proper time for an observer and the metric components scale 
as 


~g^ = 0(1)+ 0(2), (3.6) 

9oF™ = 0(1.5) and (3.7) 

4 PN = W (3-8) 

In the following, we will refer to a tensor obeying this scaling as expanded to 0(2), since 

any scalar constructed out of it will be at least 0(2). The 10 parameters 7 , (3, £, Ci and 
ai are arbitrary constants whose value depends on the specific theory in question. GR has 
7 = /? = 1 and £ = Ci = OLi = 0 . Canonical conformally-coupled scalar-tensor theories on 
the other hand typically lead to different values of 7 and /3 while the other PPN parameters 
remain at their vanishing GR values. The PPN parameters have been measured using a 
variety of different probes [37] and are all consistent with GR. One can then calculate them 
in alternate theories and use the bounds to constrain the model parameters. We will do 
precisely this below. 

4 Field solution and PPN parameters 

Here we will calculate the PPN parameters for the theory described by the Jordan frame 
metric (2.7). We will treat the scalar as a light, cosmological scalar driving the acceleration 
of the cosmic expansion. The potential V(<j>) and its derivatives are then all of order Hq, and 
we can neglect the mass of the scalar on the scales of interest. Screening mechanisms such 
as the chameleon and symmetron can render the field locally massive and thus hide it from 
local observations. In that case however, there is nothing to add to the standard conformally 
coupled screened theories in the PPN context, and we will not consider this case here. 

4.1 Preliminary Considerations 

4.1.1 Choice of Coordinates 

We begin by setting up a coordinate system in which to solve the equations. We will derive 
the field solution in the Einstein frame since the equations of motion are simplest there. The 
Einstein frame solution at zeroth-order in the CMB rest-frame is simply Minkowski space 1 
i.e. the vacuum solution of the field equations: 

d-s 2 = — dt 2 + Sij dx 1 dx 3 . (4.1) 

Note, however, that the Jordan frame metric is 

ds 2 = —N 2 dt 2 + dx 1 dx 3 , (4.2) 

where the lapse is 

N 2 = 1 - ||. (4.3) 

4 We ignore corrections coming from the FRW nature of the zeroth-order solution since these are negligible 
on solar system scales. 
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This differs from unity due to the presence of a time-dependent cosmological field cj) o- I n 
this sense, the proper time for a physical observer is not coincident with the coordinate time 
and so it is clear that even at lowest order the Jordan frame metric is not in the PPN gauge 
with this choice of coordinates. Note also that the Jordan frame energy-momentum tensor, 
2/V~9$S m /5g fll/ is related to the Einstein frame tensor introduced in the previous section 
by [16, 21]: 

V‘ v = J 1 - (4.4) 

We will consider the metric sourced by a static object of finite extent and will hence treat its 
internal structure using a fluid description. Since the Jordan-frame T^ u (rather than T^ u ) 
is the covariantly conserved energy-momentum tensor, all fluid variables such as the density 
and pressure must be defined in this frame. We will hence consider the Einstein frame as a 
calculational tool; we will not assign any physical meaning to T^, it is simply a source in 
the field equations. We will transform quantities to the Jordan-frame PPN gauge once the 
solutions at 1PN and 2PN have been found. 

4.1.2 Preferred Frame Effects 

Our final goal is to derive the PPN metric in the rest frame of the Solar System, since that 
is the frame in which observations are being made. There are two possible ways to proceed. 
First, one could directly derive the solution for the field in the Solar System frame. In this 
frame, the background field gradient, which is simply = </>o(l, 0,0,0) in the CMB frame, 
is Lorentz-boosted to 

, a . if-, w2 3 U! 4 \ 

( do4>o)ss = <j>o ( 1 + + ~8~) 

(difa o)ss = - 4>o m , (4.5) 

where we have expanded to 0(2) assuming w l = 0(0.5) following standard PPN counting. 
Once the solution is obtained, one performs a gauge transformation including terms involving 
w to obtain the metric in PPN form, Eqs. (3.1)-(3.3). One can immediately see why this 
theory predicts preferred frame effects: <j>o(t) is only isotropic in the CMB rest frame. In 
any other frame it has a spatial gradient proportional to Wi. In any scalar-tensor theory, 
an evolving background field leads to preferred-frame effects. In purely conformal theories, 
these are suppressed by powers of the ratio of the Solar System time-scale (years) to the 
Hubble time since </>o ~ Hq. These effects are usually neglected and we do the same here. 
The disformal coupling, on the other hand, adds additional explicit preferred-frame effects 
that scale as <f>o/A and thus are not suppressed if A ~ Hq. This is required to obtain novel 
effects in cosmology [23]. 

The second approach solves the field equations in the rest frame of the CMB, where the 
bulk motion of the Solar System is included in the energy momentum tensor; to lowest order, 
^CMB(^f) _ y SS^ _|_ Wi ' Q ne then performs a gauge transformation to the PPN gauge in the 
absence of bulk motion, i.e. Eqs. (3.1)-(3.3) with Wi set to zero. Finally, the PPN metric is 
boosted to the Solar System frame by a Lorentz boost (“post-Galilean transformation”), as 
described in Sec. 4.3 of [38], which reintroduces the terms involving w l in Eqs. (3.1)-(3.3). 

Both approaches are equivalent since the underlying theory is Lorentz invariant, and 
since the physical relative velocity between CMB and Solar System frames is included in 
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either case (it is merely absorbed in T^ v in the second approach). While the first approach in 
principle keeps the physics more clear, the gauge transformation from the Einstein-frame, in 
which the field solution is obtained, to PPN gauge is somewhat cumbersome in this approach, 
as many terms involving powers of w 1 and various contractions with PPN potentials need to 
be kept. The second approach on the other hands avoids these complications and so we follow 
it here. Thus, we begin by deriving the Einstein frame held solution in the CMB rest-frame, 
and then transform the resulting CMB Jordan-frame metric to the PPN form with w = 0. 


4.1.3 Ansatz for the Solution 

We expand the held as 

(t> = <Po + <t>i + <t>2, (4.6) 

where 4>o is the cosmological background held in the CMB frame, while <j>\ ~ 0(1) and 
(f >2 ~ 0(2). We expect that <fi o ~ Hq and 4>o ~ Hq in the CMB frame since the evolution of 
<f >o is driven by the cosmological background (recall that <f> is dimensionless in our convention). 
Looking at equation (2.4) with a</> = = 0 one can see that, at lowest order, there are two 

parameters multiplying the source terms for the held: 



and £ 


<t>o 

A 2 ' 


(4.7) 


We will take both of these to be small numbers and will work only to leading-order in both. 
We will see that this approximation is self-consistent once the constraints from the PPN 
parameters have been imposed. We expand the Einstein frame metric as g lw = + h )lv . 

and impose the gauge choice 


dX~\doK = -\d 0 h m and (4.8) 

<V4' - \dih^ = 0, (4.9) 


with hv = rf°h av . With this gauge choice, one can write the Einstein frame metric as ' 5 


g 00 = -l+ 2 yi + 2 x 2 , g 00 = -1 - 2xi - 2x 2 - 4xf 

90 i = B u g 0i = -B i (4.10) 

9ij = (1 + 2'L 1 )<%, r/4 = (1 - 2Ti+4T 2 ) <W, 

where our gauge choice implies 

3^ = 330^1 (4.11) 

doBk = 3^X2, ■ (4-12) 

The PPN order of the metric perturbations is xi ~ 'Ll ~ ^(1), Bi ~ 0(1.5) and X 2 ~ 0(2). 
The resulting expression for the Jordan frame metric to 0(2) is given below in Eq. (4.39). 

5 Note that the symbol \ is often used to denote a quantity referred to as the superpotential in the literature 
[38] . We will not use the superpotential in this work and use x to refer to perturbations of the 00-component 
of the metric. 
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4.1.4 The Energy-Momentum Tensor 

As remarked above, the energy-momentum tensor must be defined in the Jordan frame and 
so one has 


T^ = p 


i + - + n 

p 


u^u 1 ' + pg 


p,v 


(4.13) 


where u^ = dx^/d s, s being the proper time for an observer in the Jordan frame. The 
quantities are p ~ 0(1), the density, p ~ 0(2), the pressure and II ~ 0(1), the specific 
internal energy per unit mass, in the rest frame of the fluid. 

We will ultimately want to change coordinates so that the lapse is unity, which we can 
do by setting dT = N dt, and so the velocity measured by an observer in the Jordan frame 
is 


In this case we have 


dx l 1 dx l 

dT ~ N d t ' 

(4.14) 

dx^ _ js 

■ d ~ = (7,7 Nv), 

(4.15) 


where 7 = dt/ds and s is the proper time for an observer. In this case, one can calculate 7 
using the normalisation condition g fll/ u tJ 'u l/ = — 1 to find 

„,2 ' 


1 


N 


7=^(l + ^ + 7)+®(2). 


N 2 


Using this, we can find the Jordan-frame energy-momentum tensor to 0(2): 

2xi 

N 2 


rj-i 00 _ 


P 

N 2 


1 + n + v 2 + 2x1 


rj-tOi _ r w 

' “ N 

T ij = pv i v j + pS ij . 


pv 


(4.16) 


(4.17) 

(4.18) 

(4.19) 


We can then find the Einstein-frame energy-momentum tensor using equation (4.4). Note 
that y 7 1 - 2X/ A 2 = 1 - T/2 + 0(1) = N + 0(1). Since T 0i ~ 0(1.5) and T** ~ 0(2) one 
simply has T°* = NT °* and = NT l T T 00 contains both 0(1) and 0(2) terms. To 0(1) 
we have X = (1 — 2\'i)</>g/2 and hence 




p_ 

N 


1 + n + v 2 + 


2 xi 

N 


(4.20) 


We will also need the lowered form of the energy-momentum tensor and the trace. One finds 


Too = [1 + n + V 2 - 2xi N] , 

(4.21) 

£ 

II 

1 

4 

O 

II 

1 

(4.22) 

= T ij = N[pv i v j + p8 ij ] 

(4.23) 

T = --^[l + U + v 2 + xi T] + pv 2 N + 3 pN. 

(4.24) 


These are all the quantities that we need to compute the solution at 2PN. We will do 
this by solving the trace-reversed form of Einstein’s equation, 

RfMu = 87 tG (t iiv - ^g tw T^ , 
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and the equation of motion (2.4) for 0 to the appropriate order. Note that the contribution 
of the scalar field to the energy-momentum tensor contains terms proportional to 0o that 
are unpaired with factors of A -1 , which we neglect. Compared to the matter variables and 
terms proportional to T and S, they are suppressed by the ratio of the dynamical time of 
the system (of order a year in the Solar System) to the Hubble time H which is of order 
Id 10 years. 


4.2 Solution at 1PN 

At 1PN the only quantities we need to calculate are xi> HR and 0i- We are working in 
the Einstein frame, where due to the absence of anisotropic stress we have xi = HR. The 
00-component of (4.25) gives 




the solution of which is 


Xi = Ti = 


U 

N' 


Next, we need the field equation for (j>i, which is 

V 2 0i = 8vrGEp, 

which is solved by 

01 = -2 SC/ 

to leading-order in S. These are all of the solutions at 1PN. 


(4.26) 

(4.27) 


(4.28) 

(4.29) 


4.3 Solution at 2PN 

At 2PN we need to find the metric potentials \2 ~ 0(2) and Bi ~ 0(1.5) as well as the field 
02 - The following identities will be useful: 

U 0 i = -^-V 2 (Vi~Wi) (4.30) 

UiU't = V 2 Q[/ 2 - . (4.31) 

Note that these differ from their usual form in the literature [38] because we are working 
in a coordinate system with a non-trivial lapse. Specifically, for our Jordan-frame metric 
Eq. (4.2), the continuity equation becomes at lowest order 

N- 1 p + d i (pv i ) = 0. (4.32) 

The components Rqq and R^i of the Ricci tensor are given in many standard references (see 
[38] for example) but the calculation implicitly uses the 1PN solutions (4.27) and the identities 
(4.30) and (4.31) in the standard forms without the factors of N. The reader attempting 
to reproduce the calculations in this subsection should bear this in mind and, in particular, 
calculate these components explicitly. 

We begin with the vector Bi. The Oi-component of (4.25) is 

^V 2 Hi + = 8nGpvi. (4.33) 
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Using the solution (4.27) and the identity (4.30), one finds 


7 1 

Bi = —Vi --Wi. 

2 2 

Next, we find X 2 using the 00-component of (4.25), which gives 


(4.34) 


U 2 


1 


- ~J^2 + 2$1 + ^3 $2 + jy $ 3 + 3iY$4 ’ 


(4.35) 


Finally, we need the scalar to second order. Using the scalar’s equation of motion (2.4), we 
find 

V 2 0 2 = 4>i - E (8V 2 4>2 + 2V 2 4>3+2V 2 4>i) (4.36) 

at 0(2). The term proportional to 0 ± can be dealt with by combining the gauge conditions 
(4.11) and (4.12) to find 

3 " ”t2„ 


N U = V ^ 2 ' 


(4.37) 


Using this, one finds 




28 


(4.38) 


These are all of the 2PN quantities. 


4.4 The Jordan Frame Metric 

Now that we have all of the metric potentials and field solutions, we can calculate the Jordan 
frame metric using (4.2). We continue to work to leading-order in T and S in what follows 
and we obtain 


9oo = -N 2 1 - 


2xi 

N 2 


|2-2W) +0 {3) 


~ _ D I 0001,* , 00024 , /0/0 

goi - £>i —^ 2 —I—^ 2 —^ 0(2.5) 

9ij = (1 + 24 ' i) 8ij + 0(2). 


(4.39) 

(4.40) 

(4.41) 


This result is clearly not in PPN form: the lapse is not unity, and the last term in goo is 
0(1.5). Also, there are 0(1) and 0(2) terms in the Of-component. We need to perform a 
gauge transformation to get this into the standard PPN gauge. 

We begin by changing coordinates such that dT = N df. This leaves the 'ij-components 
unchanged but the other components are 


9tt 


9T% 


1 2 X i 2x2 20 o 0i 

+ N 2 + N 2 + A 2 

Bi 0001,j 0002, i 

~N + ^V 2_ + ^V^’ 


and 


(4.42) 

(4.43) 


thus eliminating the lapse. 








Next, we need to perform a post-Newtonian gauge transformation x^ —> x^ to bring this 
into the PPN gauge. This is a second-order transformation and so we write = x /; — ^ ^. 
where £ n is 0(n) in the PPN counting scheme. One must then expand the formula 

^(x“) = |g^A(x“(xO), (4-44) 

where a hat denotes the new metric in the new coordinate system, to second order. The 
perturbations in the new gauge can be computed explicitly using the relations given in [39- 
41]. Note that since the new metric is written in terms of the new coordinates we do not 
need to expand the metric potentials separately. 

The 0(1.5) term in goo as well as the 0(1) term in go* can be removed by choosing 

= ^- (l, o) . (4-45) 

Note that in the transformation of the metric, always induces terms of order T. Since we 
work to linear order in that parameter, any terms 0(£ 2 ) can be neglected. 

Finally, we can remove the 0(2) term in go* by performing a time shift at 0(2), 


£ = ^ {2cj> 2 - 4SG 2 ) , £ 2i = 0. 


Inserting the field solution, the metric then becomes 

~ , 2 U 2C/ 2 , 4 $! , 44>2 , 24>3 , 6$ 4 , 

50 °- _1 + iv3 _ A^ + ^ + A^ + A^ + ^ + 0{ ^> 


90i — 


1 

'N 


l v ‘ + l w ‘ 


+ 0(2.5) 

9ij = ( 1 + -jy'j ^ + 


(4.46) 


(4.47) 

(4.48) 

(4.49) 


Comparing (4.47) with (3.1), we can see that the metric is still not in PPN form, since 
the coefficient of U in 500 is not unity. The reason for this is that G = M p 2 /8n is not equal 
to the locally measured value of Newton’s constant Gn- We need to normalise the metric to 
PPN form and so we define 


Gn = w - G 


1 + -T 
2 


(4.50) 


This is the gravitational constant measured in the solar system. Next, we rescale every metric 
potential so that Gn and not G appears in their definition. That is, U is now defined such 
that V 2 t/ = —47 tGnP- Doing this, we find 


goo = -1 + 2 U- 2 N 2 U 2 + 4AT4>! + 4JV$2 + 24> 3 + 6N 2 $ 4 + 0(3) (4.51) 

71V 2 N 2 

90i = -—14; - — Wi + 0(2.5) (4.52) 

9ij = (1 + 2N 2 U) Sij + 0(2). (4.53) 

This metric is in the proper PPN form in the CMB frame, where w l = 0. Next, we proceed 
to transform it to the Solar System frame. 
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4.5 PPN metric in the Solar System Frame 

As above, we shall call the velocity of the Solar System’s centre of mass relative to the CMB 
background w l . The motion can be incorporated into our picture quite simply by performing 
a Lorentz transformation corresponding to w l on the result Eq. (4.51) in the CMB frame. 
This coordinate transformation preserves the post-Newtonian character of our metric since 
w = \JW{W 1 is the velocity with respect to the CMB dipole and hence Wi ~ 0(370 km/s) ~ 
0(0.5)). An explicit treatment of the transformation has been well-documented (Sec. 4.3 in 
[38]) and our metric takes the form 

g 00 = -l + 2U - 21V 2 1/ 2 + 41V4>i + AN<$>2 + 24*3 + 61V 2 <h 4 + 3T w 2 U + Tw i w j U ij + 4YVTA 

(4.54) 

71V 2 N 2 

90i = - Y~Vi - — Wi + T Wi U + WUij, (4.55) 

9ij = (1 + 21V 2 1/) Sij. (4.56) 

At first glance, it might seem surprising that we were able to derive the preferred-frame 
effects by working in the CMB rest frame. However, this is indeed possible because preferred- 
frame effects are present implicitly within the matter velocities v l in the CMB-frame metric. 
To see this, note that at lowest order v l can be decomposed into 

v l =^ss + wl , (4-57) 

where fg S is the matter velocity with respect to the centre of mass of the Solar System frame. 
Furthermore, notice that there are ati contributions in (and only in) the coefficients of the 
PPN potentials whose definitions involve v l . e 

4.6 The PPN Parameters 

We are finally in a position to extract the PPN parameters, which can be done by comparing 
the metric Eq. (4.54) with equations (3.1)-(3.3). We obtain the following non-vanishing PPN 
parameters: 

7 = [3 = 1 — T, «i = —4T, and 07 = — T. (4.58) 

The parameters 7 and (5 are the PPN parameters that are commonly modified in scalar- 
tensor theories. <37 and «2 parametrize preferred-frame effects due to the motion of the Solar 
System relative to the cosmological background gradient fo. Due to the relative motion of 
the Solar System with respect to the cosmological rest frame, it induces (apparent) preferred 
spatial directions locally. We will discuss some observable consequences of this below. 

The preferred-frame parameter <33 as well as the integral conservation-law parameters 
Q on the other hand all vanish, as does the Whitehead parameter £. One may worry that 
this is true to leading-order in T and £ only and that one should go to next-to-leading-order 
to derive their values. This is not the case. According to a theorem of Lee, Lightman and Ni 
[42], any diffeomorphism-invariant theory of gravity that can be derived from a Lagrangian 
is at least semi-conservative, which implies that the above mentioned parameters are zero to 
all orders. Our theory falls into this class and so there is no need to go beyond leading-order. 
The Whitehead parameter is also zero to all orders since the Whitehead potential 4*w does 

6 With the exception of the potential A. which does not appear in conservative theories without preferred- 
location effects. 
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not appear in the solutions for the metric potentials or the field. Equation (4.58) constitutes 
the main result of this section. 

In appendix A we show that the disformal contributions are unchanged when considering 
the general conformal/disformal theory A{(f>), B{cj)). Specifically, the conformal factor as = 
dinA(</>)/dln 0 is independently constrained by the light-bending parameter 7 to be a$ < 
£>(10 -3 ). This reduces any interactions of conformal and disformal effects to negligible levels. 
The parameter values given in equation (4.58) then remain valid apart from a trivial rescaling 
T — > B 2 ((j) o)T. The parameters (3^ = dlnB/dln</> and are only constrained very weakly 
(/ 8^ < G(10 5 ), < O(10 9 )) through their contributions to the PPN parameter (3 (see 

appendix A). 


5 Constraints 


5.1 A concrete model 

We have seen above that only T and not £ is constrained by solar systems experiments. 
The strongest constraint comes from 02 , which is constrained to satisfy |ck 2 | < 4 X ICE 7 
[37, 43] from limits on the Sun’s spin precession, implying |Y| < 4 x 10 - ' (we will discuss 
this constraint in Section 6 ). In order to investigate the implications for disformal gravity, 
we study a concrete representative of the disformal models considered in the cosmological 
context: 

A{<f>) = e a **, and V((j)) = mge“ A< ^. (5.1) 

Following the discussion above, we set the conformal coupling = 0 from here on. This is 
the model studied by [ 22 ] who found a dark energy dominated fixed point with 

A 3 H 2 

<f> 0 = A H 0 and =-(5.2) 

when A < \/6. Note that this implies T = —2S/A, which is consistent with our assumption 
that T ~ £. Assuming that we are close to this fixed point today, the constraint becomes 

A 2 ^ < 4 x 10- 7 . (5.3) 


The region in the A/Hq-X plane where this is satisfied is shown below in figure 1(a). One 
can see that A/Hq > 3x 10 3 is required, which gives A > ICE 30 eV. The disformal coupling 
is often stated in the form [19] 


C D 


T^d, l( l>d v (b 

M 4 


(5.4) 


in the decoupling limit. In this case, one has the relation AT 4 = Alp] 2 A? [21] and so our 
constraint translates into the bound A4 > 10 2 eV. Note that this is two orders-of-magnitude 
stronger than the previous bound using solar system tests M > 0(eV) [21], which was found 
using an estimate of the Eddington light bending parameter 7 . 

An independent constraint can be obtained from the time-variation of Newton’s constant 
in the Solar System. Taking the time derivative of equation (4.50), we find 


Gn _ 3T _ 3T 
Gn 2(1 + '^f) 2 
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Log 10 — 
Ho 


Figure 1. Left: The region in the A/Hq-X plane where the constraint from measurements of a-i is 
satisfied. 

Right: The region in the log 10 /3^,-A plane where the constraint from LLR measurements of the 
time-variation of G are satisfied. 

Lunar Laser Ranging (LLR) constrains this quantity to be less than 1.5 x 10 -13 yr~ 3 [44], 
which gives 

^AT (2/9 0 - A) < 0.002, (5.6) 

where we have used equation (5.2). This allows us to constrain as a function of A by 
simultaneously imposing the constraint (5.3). In figure 1(b) we plot the excluded region in 
the Ptf )-A plane when T just satisfies the constraint (5.3) i.e. T = 4 x 10 - '. One can see that 
is relatively unconstrained and values as large as 10 3 are not excluded. 

5.2 Implications for cosmology 

In this section we will show that the constraint > 3 x 10 3 implies that disformal effects are 
negligible in the context of the cosmological evolution. We begin by studying the expansion 
history described in the disformal model by the modified Friedmann equation in the Einstein 
frame [23] 

H = o ji r 2 ( P m Pde)i (5-7) 

oM p i 

where' 

Pm = 3 Hp m - ( ^ M ^ > ; p de = + m % e - x *\ and cf> + 2H<j) - Am^e - ^ = jf~Q ■ 

(5.8) P 

Note that we have coupled the field to all species of matter, unlike [23] who only couple cj) to 
the cold dark matter component. The disformal effects enter through the coupling Q given 

7 Comparing the definitions of the actions, (j> = f and V = M 2 ]2 Kef, where <j> re f, Kef are the corre¬ 

sponding quantities defined in [23]. Also, = B 2 /A 2 M P \ 2 in our notation. 
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by 


MpiQ 



< 10 ~ 7 pm 


(z = 0), 


(5.9) 


where we have defined (/)o such that B(cf> o) = 1 and for the order-of-magnitude estimate at 
the epoch today, we have used p m M p i 2 Hq, <j) ~ H 0 <f>, and A /H 0 >3x 10 3 . Note that the 
10“' bound was found assuming that (3$ ~ 0(1), which is the theoretically natural value. 
This parameter is only weakly constrained by PPN measurements and we have found that it 
can be as large as G(10 3 ), in which case the bound above is weakened to M p \Q < 10~ 4 p m . In 
fact, values of ftp > 0(1) lead to phantom universes and are hence strongly disfavoured [45]. 
It is then clear that the fractional modification to the background expansion Eq. (5.7) due 
to the disformal coupling is of order 1CP 7 -10 -4 or less for models that satisfy Solar System 
constraints. This is much smaller than current and near future observational uncertainties 

[ 46 ]- 

We now turn to the growth of linear density perturbations, 5 m = 6p m /p m . The growth 
equation for 5 m has been derived in [23] 8 : 


3m T 



QM P 1<A . 

/— I ^ 

V 2 Pm J 


^4vrG + 



Pm 3n 


(5.10) 


This again reduces to the GR result if Q is set to zero. We then see that the fractional 
correction to the growth factor is constrained to be less than 1CP'-10~ 4 due to our by Solar 
System constraints. In fact, the effect will be smaller numerically, as the disformal coupling 
only becomes relevant at late times, when p& e ~ p m , so that the growth is modified by order 
Q only over the last Hubble time. 

We see that for models satisfying the constraints from the Solar System derived in the 
previous section, in particular from 07 , the impact on all cosmological observables is highly 
suppressed. Finally, we point out that there are two ways to evade these constraints: first, 
one could add one of the well-known screening mechanisms to the scalar field action, i.e. the 
Chameleon or Vainshtein screening. Second, one could drop the weak equivalence principle 
and let the scalar held couple only to dark matter and not to baryons. In the latter case, the 
Solar System constraints do not apply as they were derived from baryonic objects. 


6 Conclusions 

In this work we have studied the local behaviour of disformal gravity theories to post- 
Newtonian order and have calculated the parameters appearing in the parametrised post- 
Newtonian (PPN) metric. The theory is semi-conservative and the Eddington light bending 
parameter 7 , the non-linearity parameter (3 and the preferred frame parameters 07 and all 
differ from GR, with the differences proportional to T = <^>q/A 2 . The strongest bound comes 
from a 2 , which imposes the constraint A/Hq > 3 x 10 3 . This is two orders-of-magnitude 
stronger than the previous bounds using solar system constraints [21], We were able to place 
a weak constraint on the first derivative of the disformal factor [3^ < 10 3 using Lunar Laser 
Ranging constraints on the time-variation of Gn- 

8 Note that [23] work in conformal time. We have translated their results into coordinate time. 
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The constraint on «2 comes from the near perfect alignment of the Sun’s spin axis with 
the orbital angular momenta of the planets. The term in the PPN metric proportional to «2 
leads to a torque on the Sun which induces a precession of the Sun’s spin axis, contributing 
to the misalignment of the axes of spin and planetary orbital angular momenta. Ref. [43] 
obtained a constraint on a .2 by integrating the motion of the Solar System relative to the 
cosmological reference frame over the past 5 Gyr. Caveats to this treatment given the un¬ 
certainties of galactic evolution have been pointed out in [47]. As shown there, the orbit 
of Mercury provides an independent constraint |« 2 | % 4 X 10“ 5 , with further improvements 
possible. While a detailed discussion of the «2 constraints goes beyond the scope of this 
paper, we point out that an integration time of 5 Gyr is still relatively short compared to 
the Hubble time, over which the background scalar field evolves. 

In terms of the cosmology, deviations from GR are non-negligible at the background 
and linear level when M. < OfTCP 3 eV) [23] and here we have constrained M > 10 2 eV. In 
agreement with [ 21 ], we have found that there is no “disformal screening” in these theories. 
These constraints then rule out disformal theories as a potential driving mechanism of the 
cosmic acceleration since they require a cosmology that is indistinguishable from ACDM. 
In this sense, the acceleration is due to the cosmological constant. One caveat is that the 
theory we have considered here assumes that there is a Jordan frame and so the scalar couples 
universally to all matter species. For this reason, the theory satisfies the equivalence principle. 
If one were to break this assumption and couple to dark matter only our constraints would 
be circumvented since solar system objects would not source the field. 

We end by discussing the prospects for improving our constraints using other astrophysi- 
cal probes. Binary pulsars have provided some of the most stringent tests of general relativity 
and conformal scalar-tensor theories to date [48], and so one may wonder whether the same 
is true for disformal theories. In addition to the PPN formalism, there is a parametrised 
post-Keplerian (PPK) framework for binary pulsar observations. Ref. [49] have shown that 
three of the PPK parameters can be obtained directly from the PPN parameters 7 and (5. 
These are constrained to the 10 “ 6 level at most 9 and so one cannot improve the constraints 
found here using these measurements. The most accurately measured parameter is P ppk , 
the rate of orbital decay. This is constrained to the 10 “ 12 level, however, the power emit¬ 
ted into scalar radiation typically scales as the square of the scalar charge, which scales like 
T 2 ,S 2 (see [ 21 ]) and so one expects this to yield constraints at the 10 6 level, which are not 
as strong as the bound obtained using 0 : 2 , although a more detailed calculation is required 
to confirm this. Future observations that constrain the PPK parameters to higher precision 
have the potential to improve our bounds, but, for now, they are the strongest that one can 
obtain using the properties of slow-moving astrophysical objects alone. 
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A Calculation of the PPN Parameters: General Theory 


In this appendix we will generalise our calculation in section 4 to the general theory where 
the Jordan frame metric is 


9/j.iy 


A 2 (0) 



B 2 (0) 

A 2 



(A.l) 


Now that we have gained some intuition from the previous calculation, it is possible to greatly 
simplify this without the need to repeat the entire calculation using brute force. Recall from 
the previous section that, even though we calculated the field to 2PN, neither the 1PN nor 
2PN held contributed to the PPN parameters at leading-order. The reason for this was that 
01, 2 r ^~ I E and their contribution to 500 scaled like 2 which meant they only contributed 
terms that were higher-order in S. We therefore examine the changes to the calculation that 
occur in the general case to discern whether there are any new leading-order contributions to 
the Jordan frame metric. Purely conformal contributions are already constrained by previous 
analyses of scalar-tensor theories [51] and so we are interested to see if pure disformal and 
mixed terms will yield new constraints after the calculation. 

There are two new parameters that enter in the general case: a^ = 0:^(0 0 ) and 0$ = 
0<p(4 > o), defined in equation (2.5). Expanding in these, one has, to 0(2) 

A(0) = A 2 [l + 2(f>i + (20 2 a 2 + 2(f>2a<f) + 0iQ^)] , and (A.2) 

B(4>) = 1 + 20i 0^ + 2 0 2 /^ + 2(/>20<p + 0i 0$, (A.3) 

where A 2 = A(cf> o ) 2 and we have set B(0 0 ) = 1 as before since we can absorb it into A. 
In addition to this, the definition of T and E are modified to include a factor of £> 2 (0o) so 
that T — B 2 ((j) o)T and E —> B 2 (cj) o)E. In fact, since we are setting B(cj) 0 ) = 1 this is not 
important in what follows but it does have implications for the LLR constraints derived in 
section 5 where we take time-derivatives of T. Using the same Einstein frame coordinates as 
the previous calculation, the Jordan frame metric is 


500 = ~N 2 A 2 
0001 


1 - 2 jk - 2 w - - 10 (v*<h + VIA + /?>?) 


- 2 - 


0^2002 + ^“001 + 2a^02 + 2 «00i + oi'^i — / ^a ( j ) 0 ( f ) ^^(j){ 


/ J .2 


01X1 


T 


9oi — A" Bi + A" 

9ij = A 2 [1 + 2a ( f > (f)i + 2Ti] . 


^<9*01 + ^1^*02 + 2 ( a <t> + A/O^01<9i01 


and 


(A-4) 

(A-5) 
(A- 6 ) 


If we want to find the field to 1PN we only need the energy-momentum tensor (and its trace) 
to this order. At zeroth-order, one has 7 = A -1 , where the Jordan-frame Lorentz factor 7 is 
defined below equation (4.15), and so one finds 


^00 _ _ rp _ 


A A p 

"aT’ 


(A.7) 


where we have used the general relation between the two energy-momentum tensors 


T^ v = A b (4>)\ 1 - 


2B 2 (4>)X 


A 2 


rpfljy 


(A. 8 ) 
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The scalar field equation at 1PN is then 


V 2 ^»i = 


8nGpA 4 
N 3 


f(</>o,<fro, M with f(4> 0 ,<j) 0 ,4> 0 ) = a^- T(a 0 - /^) + E, (A.9) 


so that the 1PN solution is 


A 4 f 

* = ~ 2 J^ U - 


(A.10) 


The general case differs from the minimal one in that it contains factors of T and a$ mul¬ 
tiplying U. In the purely conformal case where [G = T = E = 0 one has / = and the 
Eddington light bending parameter is [51] 


It - 1 



1 + 2 a 2 


(ATI) 


One can see from (A.4)-(A.6) that this contribution is not affected by disformal contributions 
and so, in the absence of any fine-tuning 10 or screening mechanism, the Cassini bound |y—1| < 
2.1 x 10 ~' 5 [52] imposes the constraint a$ < 10 -3 . The lapse N is unchanged in the general 
case (this is a consequence of choosing B(<j) o) = 1) and so, looking at (A.4)-(A.6) 11 , the 
only non-purely conformal leading-order corrections to the PPN parameters found in the 
minimal theory must be proportional to or ct^E 12 . Demanding that we do not fine- 

tune different contributions to the PPN parameters means that T is still constrained by the 
«2 constraint such that T, E < 10 -7 . Therefore, any additional contributions present in the 
general theory automatically satisfy the PPN constraints. In particular, there are no new 
bounds on the conformal parameter a.^ and the parameter (3$ is completely unconstrained. 
The one caveat to this is that we have assumed that Ps ~ 0(1). One can see that there is a 
very weak requirement that (3$ < 10 5 due to a contribution to 7 of 0(a ( j ) f3 ( j > T) and a similar 
constraint [3^ < 10 4 5 coming from a contribution to (3 of 0(a^/3^ 2 T). Similarly, one can see 
that there is a contribution to [3 of 0(/3^T), which imposes the weak constraint (3 ^ < 10°. 
Note that these bounds apply when a^ and T just satisfy their bounds i.e. a ^ ~ 0(1O -5 ) 
and T ~ 0( 10 -7 ). When a$ and T assume values smaller than this, and (3^ can assume 
larger values and still satisfy solar system tests. 
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